High temperature phase transition in the coupled atom-hght system in the presence 

of optical collisions 
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The problem of photonic phase transition for the system of a two-level atomic ensemble inter- 
acting with a quantized single-mode electromagnetic field in the presence of optical collisions (OC) 
is considered. We have shown that for large and negative atom-field detuning a photonic field ex- 
hibits high temperature second order phase transition to superradiant state under thermalization 
condition for coupled atom-light states. Such a transition can be connected with superfiuid (coher- 
ent) properties of photon-like low branch (LB) polaritons. We discuss the application of metallic 
cylindrical waveguide for observing predicted effects. 
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I. INTRODUCTION 



Nowadays the investigation of phase transitions in 
atomic gases represents a huge area of experimental and 
theoretical research where condensed matter and statisti- 
cal physics are closely connected to the application prob- 
lems of quantum and atom optics, e.g., in the field of 
quantum information science - see e.g. Although 
Bose-Einstein condensation (BEC) of the atoms has been 
observed in many labs, the requirement to use extremely 
low (up to microKelvins) temperatures strictly limits 
the utilization of such an effect for practical purposes. 
It provides an important reason for studying relatively 
high temperature phase transitions. Usually, such tran- 
sitions take place in coupled matter-field systems - po- 
laritons introduced many years ago for describing the 
interaction of quantized field with quantum excitations 
in the medium, see e.g. [2]. At present the evidence of 
phase transitions with LB polaritons and their superfiuid 
properties have been observed in solid statephysics with 
semiconductor microstructures - see e.g. In such 

systems polaritons can be treated as 2D gas of bosonic 
particles ~ so-called exciton-polaritons appearing in the 
sample of quantum wells inserted in a semiconductor 
(CdTe/CdMgTe or GaAs) microcavity and having effec- 
tive mass which is many orders smaller than a free mass 
of electrons. However, high (room) temperature phase 
transition for current narrow-band semiconductors seems 
to be hard to reach due to exciton ionization. At the 
same time the time of thermalization for polaritons in 
solid state structures mentioned above is short enough 
and is in the picoseconds regime Q. In this sense po- 
laritons in atomic physics could be preferable for observ- 
ing high-temperature phase transitions. Such polaritons 
potentially possess a longer coherence time. In particu- 
lar, atomic polaritons represent superposition of photon 
and polarization of two (or multi)-level atoms and can 
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be observed in various problems of atom-field interac- 
tion where long lived coherence of quantized optical field 
strongly coupled with macroscopic atomic ensemble plays 
an essential role - cf . [1, 13 ■ 

The primary step for the experimental observation of 
phase transitions and BEC for atomic polaritons is con- 
nected with achieving thermal equilibrium in the coupled 
atom-light system. The so-called optical collisions (OC) 
have been proposed recently for this purpose. The pro- 
cess of OC represents nonresonant interaction of a quan- 
tized light field with an atom in the presence of buffer gas 
particle - see e.g. Q. Although the main features of OCs 
have been investigated both in theory and in experiment 
for a long time (see e.g. 0-11.]), the thermodynamic prop- 
erties of a coupled atom-light system have been studied 
quite recently [12r,-16j. In particular, in it has been 
demonstrated for the first time that OCs lead to the ther- 
malization of coupled (dressed) atom-light states under 
the interaction of the optical field with rubidium atoms in 
an ultrahigh-pressure buffer gas cell at high (530 K) tem- 
perature. Latterly in (isj we have proposed a theoretical 
approach of a dressed-state thermalization that accounts 
the evolution of pseudo-spin Bloch vector components 
and characterizes the essential role of the spontaneous 
emission rate in the thermalization process. Although 
we have used a two-level model for describing OC process 
with buffer gas particles (see e.g. ;9.- ll]) the dependences 
theoretically obtained in [l^ are qualitatively in good 
agreement with experimentally observed results for ru- 
bidium atoms. The predicted time of thermalization was 
in a nanosecond domain at full optical power 300 mW 
and negative atom-light detuning S/2tt = — 11 THz. The 
main features of such a thermalization are connected with 
the observed asymmetry of a saturated lineshape that 
depends on the temperature of atomic gas and detuning 
S = ujl — ijjat] and uJat are frequencies of optical field 
and atomic transition respectively. 

In a number of papers [l^ [l^ authors discussed ther- 
modynamic properties of a two-level sodium atomic sys- 
tem being under OCs with buffer gas particles in another 
limit of positive atom-light detuning, i.e. for ^ > 0. In 



2 



particular, there has been shown a strong laser generation 
in the atomic ensemble interacting with cw-laser field in 
the domain of a far blue wing of spectral line under the 
frequent collision with buffer gas particles. In this sense 
both problems, i.e. lasing that occurs in atomic system 
in the presence of OC and phase transition for coupled 
atom-light states, evoke great interest and should be clar- 
ified. 

Physically the problems under discussion are very close 
to another intriguing and long stated problem of phase 
transition for a light field (or BEC of photons) - see, 
e.g., [T7l - [20j . In fact, the comparison of lasing and 
condensation phenomena for photonic, atomic and solid 
state systems as two possible ways for achieving macro- 
scopic coherence has long been discussed - see e.g. 
[23} . The first attempts to consider some analogy between 
lasing under the threshold region and second-order phase 
transition in ferromagnets have been made in pl| . Later 
a simple Dicke model that describes an ensemble of two- 
level atoms interacting with quantized electromagnetic 
field in standard laser quantum theories has been pro- 
posed for observing phase transition of photons [l7j . In 
quantum optics phase transition in such a system has 
been interpreted as a transition to some superradiant (co- 
herent) photonic state with zero chemical potential [18l |. 
At the same time it is pointed out in [Toj that in this case 
the establishment of spontaneous static field (a field with 
zero frequency) in the medium takes place. In fact such 
a phase transition creates some ferroelectric state in the 
medium. 

From our point of view there exist important circum- 
stances that must be taken into account while consider- 
ing phase transitions mentioned above. In many cases 
including usual lasers, we deal in practice with non- 
equilibrium states of the system, cf. [12, [13]. On the 
other hand, Bose-Einstein condensation phenomenon and 
phase transition under discussion have some meaning 
only in the limit of thermodynamically equilibrium state 
of a coupled atom-light system for which chemical poten- 
tial is nonzero, cf. [25|. In this sense the applicability of a 
convenient Dicke model for the phase transition problem 
should be justified in each physical case. 

Noticing, that in previous experiments [l3| the re- 
gion of high and uniform laser intensity was about 70 
fim. Physically it means that the medium is thin and the 
lifetime of photon-like polaritons discussed in this paper 
seems to be short in comparison with the thermaliza- 
tion time of a coupled atom-light system. However, the 
implementation of special metallic waveguides of various 
configurations with the length up to a few mm's in the 
OCs experiment performed recently (see [13]) allows to 
increase the time of atom- field interaction by many orders 
due to photon trapping and confinement. In fact, such 
waveguides pave the way to the investigation of the high 
temperature phase transition problem with polaritons in 
atomic optics. 

In the present paper we continue our theoretical in- 
vestigation of thermodynamic and critical properties of 



coupled atom-light states appearing due to the interac- 
tion of two-level atoms with a single-mode optical field in 
the presence of buffer gas particles. In Sec HI] we examine 
critical properties of photonic field under the thermaliza- 
tion of coupled atom-light states. In Sec lIIII we study the 
problem of a second order phase transition in the sys- 
tem under discussion. In particular, we suggest a simple 
model of polaritons describing atom-field interaction at 
equilibrium. Relying on a thermodynamical approach we 
consider a mean-field Bardeen-Cooper-Schrieffer (BCS)- 
like gap equation for order parameter - normalized am- 
plitude of photonic field. In Appendix we discuss the 
properties of a cylindrical waveguide for realizing appro- 
priate strong atom-field coupling in a single-mode regime. 
In conclusion, we summarize the results obtained. 



II. THERMODYNAMIC APPROACH FOR 
PHOTONIC FIELD UNDER THE OC PROCESS 

An optical collision (OC) is an elementary process of a 
collision between an isolated atom of sort A and a foreign 
(buffer) gas particle of sort B resulting in the emission or 
absorption of non-resonant photons - see e.g. j^- The 
inset of FigH] includes a two-level model for characteriz- 
ing OC process under consideration. Actually, in recent 
experiments [l^ - [l^ Rabi splitting frequency, atom-field 
detuning and collisional broadening are at least three and 
more orders larger than splitting frequency between hy- 
perfine levels for 5S-5P transition in rubidium il'-lines, 
cf. [2(3] . As it is demonstrated in [l^ usual pressure- 
broadened rubidium Di and D2 lines are visible at a 
moderate optical power P — 25 mW independently. In 
this limit collisional broadening can be treated for Di 
and D2 lines of rubidium atoms separately. Contrary, at 
full optical power 300 mW for which the thermalization 
of coupled atom-light states occurs the essential asymme- 
try of Di and D2 spectral lines and a common shape of 
pressure broadening curves are clearly seen in the exper- 
iment, cf. [T^ - fT^ ]. Hence, the assumption of a two- level 
approach for OC problem neglecting hyperfine structure 
of rubidium i'-lines seems to be justified, cf. [T6| . 

We suppose that an optical field with frequency ujl in- 
teracts non-resonantly with two-level atoms having fre- 
quency spacing ujat in the presence of buffer gas parti- 
cles. Physically OC means that excitation of level \b) 
is impossible due to atomic collisions only. A theoret- 
ical description of coupled atom-light state thermaliza- 
tion can be given by the density matrix approach, see 
e.g. [l3|. In general, the interaction of two- level atoms 
with a quantized optical field in the presence of collisions 
with buffer gas particles is characterized by the follow- 
ing physical parameters. First, there is Rabi frequency 
= — Srj/ilR that describes the splitting of energy 
levels in the presence of collisions with buffer gas parti- 
cles, ^R— -y/fJg -|- (5^ is a definition of the same quantity 
without collisions when = 0; fip = '2-g\/Nph is resonant 



3 



(\d l^u \ ^^'^ 
2heov' ) 

is the atom-field interaction constant, which we assume 
to be identical for all atoms; dab is the atomic dipole ma- 
trix element; V is the interaction volume; Nph is a total 
average number of photons involved in the OC process. 

The collisions with buffer gas particles are determined 
by two parameters 77 and 7 which can be expressed phe- 
nomenologically via the average phase shift that is ac- 
cumulated during the collision, cf. Q. Physically pa- 
rameters 7y and 7 can be also connected with molecu- 
lar potentials of a compound A-l-B system. Parameter 
77 depends on the difference AC/ of diagonal matrix ele- 
ments for interacting atoms A and B due to their colli- 
sions (scattering) . According to the quantum mechanical 
approach to the OC problem, AC/ can be represented as 
AC/ = C„i?~" for power law of atomic interaction; C„ 
is a constant that depends on a character of such inter- 
action; i? is a distance between the atom and perturber 
(buffer gas particle) for a given time; for more details 
see [m. 

Second, OC process is described by colhsional rate 
(collisional broadening) 7 that plays an important role 
in thermalization process. In a general case parameter 
7 is characterized by the density of buffer gas particles, 
molecular potentials for a compound system and depends 
on the value of atom-field detuning 5 - see [iH . In the ex- 
periment described in [l^] parameters 7 and 77 are of the 
order of terahertz. Skipping some details of our calcula- 
tions presented in [ISj, we are starting here with impor- 
tant results related to thermalization of coupled (dressed) 
atom-field states. 

Figure [T] demonstrates a calculated population of the 
upper state tibh that is proportional to the total intensity 
of spectral components as a function of atom-field detun- 
ing 5 taken for the interaction of a quantized field with 
two-level rubidium atoms characterized by mean reso- 
nance frequency of transition 382 THz corresponding to 
a weighted mean of rubidium i'-lines, cf. [26|. In the 
paper we focus on the so-called perturbative limit when 
atom-field detuning 5 is large enough, i.e. when inequali- 
ties \5\ > 7, 77 > r^o are held. In this limit one can neglect 
phase shift introduced by OC, assuming that 17 7?. 

At the same time the scaled atomic population of the 
upper level Ubi, in this case yields 
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where P is a spontaneous emission rate. 

First, we examine the role of atomic collisions in the 
thermalization process assuming that TS'^/jHq ^ 1; a 
blue dashed curve in FigH] corresponding to this limit. 
For negative detuning S <0 from ^ we get 
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where we also suppose that inequality 
h\d\ > keT 
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FIG. 1: (Color online) Population of the upper state ai,i, as a 
function of atom- field detuning S/2n for 500-bar argon buffer 
gas at difTerent values of the resonant Rabi frequency JIo/Stt. 
The parameters are: j/2n = 3.6 THz, F ~ 27r • 6 MHz. The 
inset includes the scheme of coUisionally aided absorption in 
a two-level atom for 5 < (left panel) and for 5 > (right 
panel) . 



is satisfied. 

The relation ^ can be recognized as a "low tempera- 
ture" limit in the framework of a convenient approach to 
existing theories for phase transitions and BEC problem, 
cf. [13]. However, in the case of OCs condition ^ can 
be achieved for high enough temperatures and for very 
large atom-field detuning S. 

The result obtained from ([2|) is in good agreement with 
theoretical predictions for the spectral line obtained for 
a far red wing (so-called adiabatic wing) - see e.g. fl]|. 
In particular, the curve of spectral intensity decays ex- 
ponentially in this case. 

Now let us switch over to a positive- valued large atom- 
field detuning, that is, to (5 > 0. In this case we get from 
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The Eq.Q reproduces a physical result for a far blue 
wing (so-called, static wing of spectral line) well known 
from the theory of OCs, cf. 111]. The result is that the 
population of the upper level is inversely proportional to 
square of atom-light detuning. 

However, according to our theoretical approach impor- 
tant features of population of the upper state are con- 
nected with the dependence of abb on temperature T of 
atomic gas. 

The thermalization of coupled atom-light (dressed) 
states occurs if condition 



p/7 < nl/s^ < 1 



(5) 



is fulfilled. In particular, for large flo the abb approaches 
Dirac-Fermi distribution function 
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In [l2|, |l3| the best result concerning coupled atom- 
light state thermalization has been experimentally ob- 
served for maximally accessible atom-field detuning 
(5/27r — —11 THz and resonant Rabi frequency 
r2o/27r = 0.1 THz (indicated by solid curve in FigH]). 
The obtained time of thermalization was 10 times shorter 
than the natural lifetime for rubidium I?-lines. 

For further processing it is important to determine a 
total atom-field excitation (polariton) density p as 
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that we consider to be constant at thermal equilibrium. 
In (O = (/V) is a normalized average photon 
number, / (/^) is annihilation (creation) operator for the 
photons absorbed (or emitted) due to atomic collisions at 
the equilibrium, A'^ is a number of atoms. 

The thermodynamic property of a coupled atom-light 
system depends on the contribution of photonic and 
atomic parts in Eq.(l7]) which is usually considered in 
some specific limits. 

The so-called low density limit p ^ 0.5 for atom-field 
excitations (polaritons) implies that 



A < 1, CTfeb < Ga 
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which is obtained at negative atom-field detuning 5 < 
at "low temperature" limit ([3]), when atoms mostly 
populate their ground state \a). 

Relation (|8^) represents the necessary prerequisite for 
observing a second order phase transition in vario us p hys- 
ical systems being at low temperatures, cf. 0, [H^ [27| . 
As for the problem of OCs a low density limit (|5]) can 
be achieved at high temperatures but for a very large 
atom-field detuning \d\. An experimentally accessible rel- 
ative photon number involved in the atom-field interac- 
tion and estimated for Rabi frequency J7o/27r = 0.1 THz 
is Nph/N ~ 6.24 X lO'^ - see [ll, [l3|. We assume that 
inequality ([S^) is fulfilled for the system under discussion. 

The saturation of atomic population is achieved at ex- 
citation density p « 0.5 for H\S\ ^ ksT. This is the 
so-called secular approximation (see Q) for which OCs 
tend to equalize dressed-state populations which is not 
the case of our study here. 

For positive detuning {S > 0) from ([7]) and ([8^) we 
obtain p > 0.5 {abb > 0.5) that corresponds to the limit 
of population inversion in a two-level atomic ensemble, 
cf. [i^- Physically, such a system behaves unstable be- 
cause of spontaneous emission processes from the upper 
level. In particular, relevant population abb diminishes 
at a large positive detuning for any fixed laser intensity, 
see FigHJ 

To find photonic field properties under OC process at 
equilibrium we consider ([7]) as an equation for order pa- 
rameter A at the given density p. Putting in ([7]) A = 
for critical value ac of vital parameter a = hS/kBT that 
determines the phase boundary between normal and "su- 
perradiant" states one can obtain: 




FIG. 2: (Color online) Order parameter A versus vital pa- 
rameter a (normalized atom-light detuning 5) for T — 530 
K temperature of atomic gas. In the inset a phase bound- 
ary for critical parameter ac versus LB polariton density p is 
presented (see the text). 



Furthermore, one can get from (j6]) - ([8]) an expression for 
order parameter A(q;) as a function of atom-field detuning 
S or temperature T 



X{a) = Ao 



1 



l + (l/p-l) 



a/ac 



1/2 



(10) 



ac = -In [(1 - p) Ip\ 



(9) 



where Aoo = \fp is an order parameter at "zero temper- 
ature" limit (d]). 

Eq.(|10|) gives an opportunity to interpret critical prop- 
erties of the photonic field occurring due to OCs as a 
result of atom-field thermalization. To be more specific 
we suppose atom-light detuning 5 to be negative. While 
modulus |5| is small enough (such as h\6\ <^ ksT, 7), 
the elementary processes of emission (or absorption) of 
photons by atoms happen independently on the atomic 
collisions with buffer gas particles during their free mo- 
tion. The atomic collisions lead to the dephasing of emit- 
ted radiation. In this sense they are worse and we deal 
with a "normal" (incoherent) state for photons. 

A physical picture is changed significantly by increas- 
ing atom-field detuning S at high buffer gas pressures. 
For large \S\ under the condition of ?i|(5| > fcsT it is no 
longer possible to ignore correlations between elementary 
acts of atomic collisions with buffer gas particles and pho- 
ton emission (or absorption). This is a case of OCs. Each 
collision happens in a very short time period comparing 
with the time interval separating two collisions. Watch- 
ing OCs for a long time when many collisions happen 
and frequent transitions of atomic population between 
dressed states occur it is possible to create some pop- 
ulation of coupled atom-light states in thermodynamic 
equilibrium due to the thermalization process. This pop- 
ulation can be established for "usual" (bare) atomic levels 
only for large values of atom-field detuning. In this case 
Eq.(|8l) represents transition to some ordering ("superra- 
diant") state as a result of forming macroscopic "spon- 
taneous" polarization for a coupled atom-light system. 
It should be noticed, that the existence of some small 
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residual polarization of the atomic system has been pre- 
dicted in 13] as a result of the atom-light (dressed) states 
thermalization. 

Figure[2]demonstrates phase transition behavior for or- 
der parameter A for the system under discussion. From 
practical point of view it is much easier to vary atom- 
field detuning S instead of atomic gas temperature T. 
The dashed curve corresponds to density p = 0.27 in ac- 
cordance with experimentally achieved detuning S/2'k = 
— 11 THz. The solid curve is relevant to low density limit 
([8]). The critical value \ac\ of parameter a is increased by 
decreasing excitation density p - see inset in Figl^l In the 
low density limit ([8]) Eg. lfTO)) is simplified and looks like 

A (a) ~ Aoo 1 — (p)" . Phase transition occurs 

for large enough detuning (or at lower temperatures). 

For positive detuning {5 > 0) the problem of phase 
transition under discussion is sophisticated. Formally, 
full thermalization of coupled atom-light states and rele- 
vant phase transition can be achieved at infinite resonant 
Rabi splitting frequency fio and atom-field detuning S, 
when the population of the upper level is a^b = 1 - see 
However, an approach to OCs that characterizes 
thermalization process is valid only for limited values of 
detuning such as \S\ <^ Uat- For any finite f2o the varia- 
tion of atom-field detuning 5 or gas temperature T drives 
a coupled atom-light system out of thermal equilibrium 
even if such an equilibrium (or quasi-equilibrium) has 
been initially achieved - Fig[TJ In this sense, the ful- 
fillment of condition ([S]) for coupled atom-light states 
permits to create a population inversion in a two-level 
atomic ensemble and to realize non-equilibrium (or quasi- 
equilibrium) transition to lasing as a result, cf. [l5l[l6l|. 



III. PHASE TRANSITION IN A COUPLED 
ATOM-LIGHT SYSTEM UNDER EQUILIBRIUM 

While OCs result in thermalization of coupled 
(dressed) atom-hght states, the Eq. lfTOl) describes the 
phase transition under discussion only qualitatively. 
Strictly speaking, a rigorous thermodynamical approach 
to critical properties of a coupled atom-light system at 
equilibrium is needed. 

We describe the interaction of a single-mode light field 
/ with an atomic ensemble by relevant (Dicke) Hamilto- 
nian - cf. [13, El- 

(11) 

where k = g^/N is a collective parameter of atom-field 
interaction, S-^j is a transition operator for j-th atom, 
Sz,j is an operator of atomic population imbalance. 

From the mathematical point of view for the definition 
of S-,j and Szj operators it is possible to use Pauli spin 
matrices representation (cf. [l8|) or explore annihilation 
(flj, bj) and creation (aj, 6^) operators for bosonic atoms 



ajbj and Sz.j 
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at the ground \a) and excited |6) states in the second 
quantization representation respectively. The latter one 
enables definitions 
see [H. 

Since the atomic medium is very dense (uat — 
lO^^cm^^) one can neglect inhomogeneous (Doppler) 
broadening because Doppler broadening doesn't restrict 
the thermalization process and is by two orders smaller 
than parameter k, that is k/27t « 0.624 THz for resonant 
Rabi frequency no/2n = 0.1 THz - cf. [11,113. Hence, 
the so-called strong atom-field coupling condition that 
plays an important role in the problem of phase transi- 
tion with polaritons is fulfilled in our case, cf. [i-Q. In 
this limit one can introduce collective atomic ladder op- 
erator S- and operator of atomic population imbalance 
Sz as follows: 



N N 

S- = E ^-,3^ — E "^^j- 



Operators introduced in obey SU(2) algebra com- 
mutation relations 



S , S 



[Sz,S_] 



'2S-. 



(II3i,b) 



For describing excitations of a two-level atomic sys- 
tem we define annihilation {((>) and creation (0^) excita- 
tion operators by u sing the so-called Holstein-PrimakofF 
transformation (cf. |29[) 



_ = y^N - 0t00, si = (j)WN Sz = 2^U-N. 

(IIli,b,c) 

The relations (|13p are preserved when operators (j) and 
^ obey usual commutation relation for the Bose-systcm: 



= 1. 



(15) 



For a large number of atoms N it is possible to treat 



operators S- and Sl_ as 



S- = Vn<p - 



2VN 



, Sl = VN(P^ - 



2VN 



The last term in (1161) characterizes nonlinear effects 
for atomic excitations under the atom-field interaction. 
However, in the low density limit defined as 



(17) 



we can obtain from ()16p simple expressions for operators 
(f) and (pi i.e. we get 



P 



1 



5_, 



At 



N 



Taking into account Eq. p4t :). it is easy to see that 
condition pT|) exactly corresponds to the relation ([SJd) 
obtained for large and negative atom-field detuning S < 
0. In this case collective atomic excitations described by 
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operator (f> characterize a macroscopic polarization P of 
atomic system. 

Taking into account Eqs. ([T4h '). (|17p . Hamiltonian (fTTj) 
can be easily reduced to the form 



(19) 



At the steady-state for a coupled atom-light system 
Hamiltonian ([T9|) can be diagonalized with the help of 
unitary transformations 



where 2 = ^ ^1 ± 



^ s'^+iK^ I Hopfield coefficients 
that obey normalization condition ^1 + '1)2 = I- 

The annihilation operators 2 in Eg.d^Dk.b) charac- 
terize two types of quasiparticles due to the atom-field 
interaction, i.e. upper and lower branch polaritons. By 
using relation (ITSI) it is easy to show that operators de- 
fined in (|20b ..b) obey commutation relations 



^„<i>\\ = k,<i>il =1 



(21) 



Below we focus on photon-like polaritons necessary for 
explaning photonic phase transition for the light field in 
the system under discussion. 

The equation for order parameter A at equilibrium 
can be derived from a variational (thermodynamic) ap- 



proach, see, for example, 18|. In their calculations the 



authors use a canonical ensemble with chemical potential 
II — applied to Dicke Hamiltonian (|lip . However, this 
approach cannot be justified even for pure photonic sys- 
tem posing phase transition at equilibrium. For example, 
as it is shown in [s^ the chemical potential for photonic 
gas being under BEC condition is nonzero. 

In the paper we use a grand canonical ensemble with 
finite chemical potential and coherent basis for the pho- 
tonic field for calculating partition function Z{N, T) — 

Tr (^e"-^'/'=B^^ ; H' ^ H - ^N^x is a modified Hamifio- 

N 

nian; Nf.x = /^/+^ i^ ^^'^ number of excitations. 

Evaluating partition function Z{N, T) under the mean 
field approximation one can obtain (cf. [1^): 

^2^tanh(^ (^2^ + 4.^2^2) 
'^lX = : —-T7^ , (22) 



,1/2 



where we introduced denotations uj^ 

Equation (1221) is similar to the gap equation that char- 
acterizes a second-order phase transition for various sys- 
tems in solid-state physics, see e.g. For finding 
chemical potential /i in our case one can use polariton 
number operator Npoi — + $2*^2- With the help of 
Eas. (fT2b ). (fMb ) and (l20l) Npoi can be represented as 



pol 



N 



N 



Nex. 



(23) 



Density p of total excitations occurring in a closed atom- 
light field system defined in ([7]) can be interpreted as 
a polariton number density, i.e. p = Npoi/N. In the 
thermodynamic limit it approaches 



oj^f tanh 



(2feBT i^at 



A2 



1/2 



1/2 



(24) 

The last term in the brackets of Eq. (fM)) characterizes 
normalized population imbalance = S^/N at thermal 
equilibrium. In particular, it is (cf. (25j): 



"^at tanh ^ 2k^ 



(2/cBT (' 



4^2 A2 



1/2 



4k2A2)^/^ 



(25) 



Combining and for chemical potential /x we 



get 



Ml, 2 = 2 i^at + ± flR^eff] , 



(26) 

where flu^eff ~ \J ^"^ ^ (p — A2 — i) is a new effec- 
tive Rabi splitting frequency. 

At low polariton densities (see ([8])) Eq. ([M)) defines 
normal state (A = 0) for upper {pi) and lower {^2) polari- 
ton branch frequencies respectively. Inserting (1261) into 
(I22]) we arrive at 



' tanh 



■^Ll,2 



(2^(-^*l,2+4.^A2)^/^) 



'^atl,2 



4k2A2 



,1/2 



(27) 



that are BCS-like equations for upper (index "1") and 
lower (index "2") branch polaritons respectively; w^^i 2 — 

i [(5 =F ^R,e}f\^ ^atia = \[-5T ^R.eff]- 

The critical temperature Tc of phase transition (for 
given atom-field detuning 5) can be obtained from (j27p 
for A = and looks like 



Tci,2 — 



atl,2| 



2kB tanh"' [± (2p - 1)] 



(28) 



The results, obtained from Eq. ([28|) are consistent with 
existing theories of BCS-like phase transition in polariton 
system, cf. [7l.[25l|. 

For very large atom-light detuning, i.e. for \S\ ^ k, a 
normalized population imbalance (|25p looks like 



Si""^ ~ -^tanh(S|5| /2kBT). 



(29) 



Equation ([29]) corresponds to equilibrium population of 
the upper level achieved at full thermalization of coupled 
atom-light states in the presence of OCs and described 
by Eq.(l6]). In the same limit the critical temperature Tc 
approaches 



Tr 



hS 



2fcBtanh"'(2p- 1) 



(30) 
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Ea. (l30l) immediately comes from Eq.0 taken for order 
parameter A = 0. 

At the same time Eq. (f28|) defines critical value ac of 
a-parameter (normalized atom-field detuning S) for fixed 
temperature T of atomic gas for which we get 



ac 



In 



1 



klTHn{{l-p)lpy 



(31) 



Equation (j31[) characterizes phase boundary between 
photonic superradiant (coherent) and normal states 
taken at A = for the given temperature T of an atomic 
ensemble. In accordance with experimental conditions of 
p^ - [T^ . fcsT ^ Hk] and the last term in (|31l) vanishes. 
In this case we recognize Eq.Q for critical parameter ac 
that characterizes OC process at full equilibrium. 

For negative atom-light detuning, i.e. for 5 < 0, we 
can interpret the phase transition under consideration as 
a transition to superfluid state of LB photon-like polari- 
tons ($2 — —/) occurring under OC process in the low 
density limit In this case order parameter A can be 
recognized as a polariton wave function. In Appendix we 
represent the main features of cylindrical waveguides for 
observing phase transition with polaritons. Superfluid 
properties of polaritons can be elucidated taking into ac- 
count a weak interaction between them and can be de- 
scribed in the momentum representation, cf. Q • Such an 
interaction arises due to nonlinear effects in atom-field in- 
teraction, described by second terms in expressions (jl6p 
for 5_ and 5*1, cf. [Slj. We expect, that in the future 
it will be possible to examine superfluid properties of 
such polaritons by using very thin-walled (up to 60 nm) 
waveguides for the experiment with Josephson junctions 
at high enough temperatures. 



behavior of photon-like LB polaritons under the low den- 
sity limit ([5]). In this respect, we discuss special metallic 
waveguides for providing necessary atom-field interaction 
and achieving suitable polariton lifetime, as a result. In 
particular, phase transition under discussion can be con- 
nected with transition to superfluid (coherent) state for 
such polaritons characterizing superposition of a quan- 
tized optical field and a macroscopic atomic polarization. 
In the presence of phase transition the atomic polariza- 
tion evolves in time with the frequency that is equal to 
the chemical potential p2 ~ ^at — \5\ for LB polaritons. 
This property enables to observe the phase transition un- 
der discussion experimentally. 

We note that the results obtained above can be useful 
for observing true BEC occurring with atomic polaritons. 
In this case we need to realize certain trapping potential 
for LB polaritons. One simple way for that is to use a 
biconical waveguide cavity for which waveguide radius R 
varies smoothly along z-coordinate. A similar dielectric 
(bottle like) cavity based on tapered optical fiber has 
been recently proposed for improving atom-light coupling 
strength, see e.g. [s^]. In our case such a cavity enables 
to create an appropriate trapping potential for photons 
as well as for polaritons in z-dimension. It is important 
that in this case the lifetime of photon-like polaritons can 
be long enough and is determined by the cavity Q-factor. 
These problems will be a subject of intensive study both 
in theory and experiment in forthcoming papers. 
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91350 and by Russian Ministry of Education and Science 
under Contracts No. n569, 14.740.11.0700. A. P. Alod- 
jants is grateful to M. Weitz and F. Vewinger for hospi- 
tality and valuable discussions. We are grateful to the 
referee for his comments. 



IV. CONCLUSIONS AND OUTLOOK 

We consider a phase transition problem in a two- 
level atomic ensemble strongly interacting with an op- 
tical field. The main features of such a phase transition 
are connected with a coupled atom-light state thermal- 
ization occurring due to the OCs process with buffer gas 
particles. Such a thermalization can be achieved exper- 
imentally for a large value of negative atom-field detun- 
ing 5. Using a thermodynamic approach (partition func- 
tion) we established a gap equation under a mean field 
approximation for order parameter A - normalized aver- 
age optical field amplitude. In the paper we represent 
simple arguments confirming that the obtained thermal- 
ization of coupled atom-light states leads to the photonic 
phase transition to the superradiant (coherent) state and 
is characterized by some ordering (equilibrium state) for 
two- level atomic system. Nontrivial solution of Eg. ([2^ 
with A ^ leads to the appearance of a macroscopic sta- 
tionary polarization of the atomic medium that is pro- 
portional to the order parameter A. 

The equilibrium properties of a coupled atom-light sys- 
tem offer a promising approach to studying the critical 



Appendix A: PHOTON CONFINEMENT IN 
CYLINDRICAL WAVEGUIDE 

Let us briefly discuss the possibility to observe the 
phase transition under discussion with LB polaritons 
in metallic cylindrical waveguide. The implementation 
of metallic microtubes for investigating coupled atom- 
light state thermalization under OC has been demon- 
strated in [ij]. The properties of an optical field in 
the ideal empty waveguide without any losses are well 
known - see e.g. (ssj . In particular, wave vector com- 
ponent fcj_^,„p = grap/R of the field in the cross section 
of the waveguide is orthogonal to z-axis and quantized; 
9mp represents the root of the equation Jmik±R) = 0, 
i? is a waveguide radius, m and p are integer numbers 
characterizing azimuthal and transversal field distribu- 
tion respectively. At the same time there exist continuum 
modes with wave vectors kz in z direction. Physically it 
means that it is possible to establish a dispersion relation 
for photonic field in the waveguide as 

UJL ^ cfc_L,mp + hkl/2'mphj (Al) 
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where we have introduced photon mass rUph = hk±^rnp/c. 
Thus, for a cyhndrical waveguide spatial degrees of free- 
dom in X and y directions are suppressed and photons 
remain confined in the plane which is perpendicular to 
2;-axis. At the same time Eq. ljAll) implies the exis- 
tence of a finite polariton mass in the waveguide. For 
example, mass nipoi of photon-like polaritons is about 
2.6 X 10~^^ kg that implies a high temperature phase 
transition in our case, cf. 

It is possible to define from Eq. (|Aip the so-called cutoff 

wave number kmp = k±^mp = gmp/R for which a waveg- 
uide mode propagation constant vanishes. In the exper- 
iment it is preferable to use a fundamental TMoi-mode 
with quantum numbers m = 0, p = 1 {gm — 2.4048) as 



a cavity mode for the atom-field interaction. In partic- 
ular, we require the fulfillment of conditions /xi_2,wl > 

ckg^' for characteristic frequencies that describe a cou- 
pled atom-light system. The condition under discussion 
can be represented in some other way by using radius 
R of the waveguide as ^°-\s\ < R < ^zrT+\5\^ where 
gii ~ 3.8317. The fulfillment of the last inequality guar- 
antees that only a single waveguide mode effectively in- 
teracts with the atomic ensemble at a large atom-light 
detuning S. Notably, that according to the condition 
presented above the diameter of the waveguide should be 
approximately of the order of wavelength A^h ~ 785 nm 
for rubidium atom _D-line transition, cf. [l^ . 
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